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A new parameterization for the dark energy equation of state(EoS) is proposed and some of its cosmological 
consequences are also investigated. This new parameterization is the modification of Efstathiou' dark energy 
EoS parameterization. w(z) is a well behaved function for z » 1 and has same behavior in z at low red- 
shifts with Efstathiou' parameterization. In this parameterization there are two free parameter wq and vv„. We 
discuss the constraints on this model's parameters from current observational data. The best fit values of the 
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I. INTRODUCTION 

In recent years, the discovery of accelerating expansion 
of the universe is an amazing development. It was firstly 
discovered by observing type la supernova (SNe la) 111 |2|], 
which can be used as standard candles L3„,4J. The cosmic mi- 
crowave background (CMB) measurements from Wilkinson 
Microwave Anisotropy Probe (WMAP) p] and the large scale 
structure survey by Sloan Digital Sky Survey (SDSS) |0 ^ 
confirm this accelerating expansion universe model. There 
are two kinds of ideas, i.e. the existence of the dark energy 
or modifications of the gravity theory, to explain this concept. 
The first scheme is most popularly discussed, and many mod- 
els have been proposed, such as the holographic dark energy 
models |8] and the Chaplygin Gas |9] . In addition there are 
also many modified gravity models, such as the brane world 
flQl and f{R) UUJ and so on. 

The dark energy equation of state(EoS) which is the ratio 
of pressure to energy density, is a prefect quantity to study the 
behavior of dark energy. If EoS is a constant - 1 , dark energy 
is the A Cold Dark Matter Model(LCDM), and maybe dark 
energy is vacuum energy. Otherwise the dark energy would 
be dynamical scalar field, such as the Quintessence I.1Z1 . the 
quintom lllSll . 

There are several way to explore the behavior of dark en- 
ergy EoS. The most popular way is to build a functional form 
for EoS in terms of some free parameters. Lots of EoS param- 
eterizations have already been discussed in the literature (such 
as 1,14.-24.1 and Refs. therein). Another way is picking a sim- 
ple local basis representation for w{z) (bins, wavelets), and es- 
timate the associated coefficients [25-27]. In addition, there 
are also some nonparametric way |28, 29|. 

In this paper, we consider a new parameterization for the 
dark energy EoS. In[20], the author developed a new dark en- 
ergy EoS parameterization: w{z) - wo + w„ ln(l + z). It is ob- 
vious that when ^ — > oo, wiz) has poor behavior and becomes 
infinite . This dark energy EoS parameterization can only de- 
scribe the behavior of dark energy when z is not very large. 
To avoid this problem, we consider a new dark energy EoS 
parameterization, which is the modification of Efstathiou' pa- 



rameterization: w(z) - wq + Wc, ln(l + ^). The value of w{z) is 
Wo at present and w{z) becomes to (wo + Wa* ln2) when ^ ^ oo. 
In this model there are three parameters in all, which is wo, Wa 
and 0,„. As shown in[30], this EoS will get to a nonphysical 
value in the far future time when redshift z approaches -1, 
namely, |w(z)| will grow rapidly and diverge. 

In this paper, we perform a global data fitting analysis on 
this new dark energy EoS parametrizations, and present con- 
straints on the model parameters from the current observa- 
tional data, including the seven-year WMAP data, Baryon 
Acoustic Oscillations (BAO) data. Observational Hubble data 
and SN Union2 sample. Since dark energy parameters are 
tightly correlated to some other cosmological parameters, 
such as the matter density parameter Q,„ and the Hubble con- 
stant //(), it is necessary to consider a global fit procedure in 
the investigation of the dynamical dark energy. The paper is 
organized as follows: In section II, we review the new dark en- 
ergy EoS parameterization. In section III, we describe Current 
Observational Data we used. In section IV, we perform the 
cosmic observation constraint, the results are also presented. 
The last section is the conclusion. 



II. NEW PARAMETERIZATION 

Let us start by presenting some of the most investigated EoS 
parameterizations(see also 1I22I1 for other parameterizations): 

' wo + WaZ (redshift) LiS-iT] 

Wo + WaZ/(l + z) (scale factor) [18. 191 



w(z) 



Wo + Waln(l -Hz) (logarithmic) [201 



where wo is the current value of the EoS, and w^ indicate the 
revolution of the dark energy EoS. The value of these parame- 
ter is determined by the observational data. When Wa -0 and 
Wo = -1, the dark energy model becomes to LCDM model. 

The first parameterization represents a good fit for low red- 
shifts, but has serious problems to explain high-z observations 
since it blows up as exp(3WflZ) when z > 1 and Wa > 0. For 




FIG. 1; the CMB anisotropy spectrum. Red line: wq = -0.7, w„ = 
-0.2; Green line: wq = -1.2, w^ = 0.2; Black line: wq = -1, Wn = 0. 



example, it can not explain the estimated ages of high-z ob- 
jects |31]. The second one solves this problem, since w(z) 
is a well behaved function for z » 1 and recovers the lin- 
ear behavior in z at low redshifts. The latter was introduced 
by Efstathiou 1.20.1 . It was built empirically to adjust some 
quintessence models at z < 4. When z approachs infinity, w(z) 
has poor behavior and becomes infinite and this is unnatural. 
Similar to the second model, let us consider the following EoS 
parameterization 



w(z) - Wo + Waln(l + 



l+z 



(1) 



•w„ 



where the value of w(z) is Wq at present and w(z) — » (wo - 
ln2) whenz — > oo. 

If there is no interaction between dark energy and other 
component of the universe, one can show from the energy con- 
servation law [p - -3a(p + p)/a] that the dark energy density 
evolves as 

Pdb(z) = Pad - Qm) exp ^ f i'^ + ^(^O] y^ j , (2) 

where pc is the critical density, it is defined by the following 
equation 



Pc = 



3Hl_ 
SjtGn 



(3) 



In this new parameterization, it is hard to write out the analy- 
sis formula of this quantity. It is calculated through numerical 
method. 

In order to study the evolution of cosmological perturba- 
tions, we use the public Parameterized Post-Friedmann (PPF) 
package developed by Wayne Hu(see e.g. |32] for detailed 
discussions on PPF method) to calculate the CMB anisotropy 
spectrum. In Figure [1] we plot the anisotropy spectram for 
different choices of wq and Wa- We observe some obvious dif- 
ferences with respect to the LCDM case on larger scales (mul- 
tipole number / < 10). And for the model wq - -0.7 and Wg - 
-0.2, there is a slit deviation on the peaks of anisotropy spec- 
trum. 



III. CURRENT OBSERVATIONAL DATA 

In order to test this new model, we use the most recent ob- 
servational data currently available. In this section, we de- 
scribe how we use these data. 

A. Type la Supernovae constraints 



We use the 557 SNe la Union2 dataset fSS"]. Following | 
13511 . one can obtain the corresponding constraints by fitting 
the distance modulus fi(z) as 



jj.,h(z) = 5 logio[£)t(z)] + no 



(4) 



where jiq = 42.38 - 5 logjg h, and h is the Hubble constant Hq 
in units of 100 km s ' Mpc ', In flat universe, the Hubble-free 
luminosity distance Di - HqcIl is 



Dl{z)^(1+z) 



I 



dz' 
Eiz'Y 



(5) 



where £(z) = mz)/Ho. 

For the SN la dataset, the best fit values of the parameters 
can be determined by a likelihood analysis, based on the cal- 
culation of 



2 
A'SN 



557 



[//thfc)-j»obsfc)]^ 

crHzi) 



(6) 



B. Baryon Acoustic Oscillation constraints 

The BAO data come from SDSS DRV |36]. The datapoints 
we use are 



rs(zd) 
Dv(0.275) 



= 0.1390 ±0.0037 



and 



Dv(035) 
Dv(0.2) 



1.736 ±0.065 



(7) 



(8) 



where the spherical average gives us the following effective 
distance measure 1.37,1 . 



Dv(z) 



p" dx 
Jo W) 



)l H(z) 



1/3 



(9) 



and rs(zd) is the comoving sound horizon at the baryon drag 
epoch. Also, Zd can be obtained by using a fitting formula 



Zd = 



1291(Q„,/z^)' 



2\0.251 



2\b2 



1 + 0.659(Q,„/i2)0 



828 



[\ + b,(^bhy'], (10) 



with 



bi = 0.313(Q,„/j2)-°-4''^[l + 0.607(Q„/zy'^4], (11) 
b2 = 0.238(Q,„/!2)°-223. (12) 



The function r^iz) is the comoving sound horizon size 
^'/<'+^) da 



rsiz) 



-I 

V3 Jo 



a^H(a) yjl + (3Qfc/(4n» 



(13) 



where Qy = 2.469 x IQ-^h-^ for Tcmb = 2.7257*:. 
So the x^ for the BAO data is given by 



2 ^ / rvfa)/£'v(z^ 0.275)- 0.1 390 r 



0.0037 

Dv(z^0.35)/Dv(z^0.2)-1.736 \^ 

0.065 I ■ ^ ^ 



C. Cosmic Microwave Background constraints 



The CMB shift parameter R is provided by 135 

™ dz' 



R(Zrec) = — ^^sinn[ Vl^ 



H(z') 



(15) 



where sinn(;i(;) is sinh(x) for Q^ > , jc for Qj. = 0, and sin(x) 
for Q^ < 0, respectively. Here, the redshift Zrec (the decou- 
pling epoch of photons) is obtained by using the fitting func- 
tion 1,4^ 

Zrec = 1048 [l + 0.00124(D.bh^r'^™] [l + gi(£ljt^y'-] , 

where 

gi = 0.0783(Qi/!2)-°-23**(lH-39.5(Oi/!T^''^) ', 
g2 = 0.560(l+21.1(Q/,/i2)i-8i)"'. 

In addition, the acoustic scale is related to the distance ratio 
and is expressed as 

Ia = -^-^sinn[ ,^ r -g-]. (16) 



FoUowing Ref. MM, the;^'^ for the CMB data is 

xImb = (^f - xf')(C-%(xf - xf), (17) 

where x,- = {lA,R,Zrec) is a vector and (C"'),j is the in- 
verse covariance matrix. The seven-year WMAP observations 
MTIl give the maximum likelihood values: lAiz^c) - 302.09, 
R(Zrec) = 1-725 and Zrec = 1091.3. In Ref. liill, the inverse 
covariance matrix is also given as follows 



ic-') 



( 2.305 29.698 -1.333 ^ 
29.698 6825.270 -113.180 
-1.333 -113.180 3.414 



(18) 



D. Observational Hubble data (OHD) 



The Hubble parameter can be writen as the following form: 



H{z) = - 



1 dz 



(19) 



\+zdt 

So, through measuring dt/dz, we can obtain H(z). In [42] and 
143. 441 . the author discuss the possibility of using absolute 
ages of galaxies to calculate the value of dt/dz- In ll44ll . the 
galaxy spectral data come from the Gemini Deep Deep Sur- 
vey 0451] and archival data [46-51]. Detailed calculations of 
dt/dz can be found in MM , and we do not discuss them here. 
Currently, we have a set of 12 values of the Hubble parameter 
versus redshift in total(see table 2 of 11521] ). The data calcu- 
lated by this way are less sensitive to systematic errors which 
is a great advantage ]|53l]. 

We can use these data to constraint different kinds of dark 
energy models and modified gravity models by minimizing 
the quantity 



2 
/VOHD 



12 

Z 



(H(zj) - H,Uzj)y 



(20) 



This test has already been used to constrain several cosmolog- 
ical models llUJel. 



IV. RESULTS 

In our analysis, we perform a global fitting to determine 
the cosmological parameters using the MCMC method. In 
our calculations, we have taken the total likelihood function 
L oc e^^'^^ to be the products of the separate likelihoods of 
SNe la, BAO, CMB and OHD. Then we get;^'^ ^^ 



X 



2,2 ,2 ,2 

■XsN '^ Xbao '^ Xcmb '^ Xohd^ 



(21) 



The results on the best fit values of the cosmological 
parameters with \cr confidence-level regions are: Q,„ = 
0.2735!»:»[^i, wo = -1.0537!«:[«2 and w, = 0.2738!°^. 
The nuisance parameters Ho used in the analysis is actually 
not model parameters with significant meanings, so we do not 
list it. 

In Figures |2] we also show the parametric spaces wq - w^ 
and Q„, - wo that arise from the joint analysis described above. 
We note that the result is consistent with the LCDM(wo = - 1 
and Wfl = 0) model in the Icr CL. To acquire more information 
on the property of dark energy, in Fig. |3]we plot the evolution 
of the EOS w(z) along with z- The following is some discus- 
sion of these findings: 

1. The best-fit results are: Q,„ = 0.2735, wo = -1.0537 
and Wfl - 0.2738. Note that here the results are max- 
imum likelihood values. The value of w(z) is -1.0537 
at present and w(z) equals to -0.8639 when z ^> o°. 
We find that the best-fit dark energy model is a quintom 
model [O], whose w(z) crosses the cosmological con- 
stant boundary w = - 1 during the evolution. And the 
redshift is 0.2766 when w(z) crosses the cosmological 
constant boundary w = - 1 . 




is in the band {-1.7, -0.4), which is relatively narrow. 



-1.6 -1.4 



-0.8 -0.6 




4. The best fit value and Icr confidence-level regions of 
Q,„ is 0.2735;^yy[g3, which is also consistent with the 
constraint on Q„, in the LCDM model and the CPL 
model. 





FIG. 3: The evolution of w(z) along with z for the model considered 
in this manuscript. The result is consistent with the cosmological 
constant in the Icr CL. 



V. CONCLUSION 



FIG. 2: Marginalized probability contours at Icr and 2cr CL in 
the n„, - Wo and wq - w„ planes for the model considered in this 
manuscript. The result is consistent with the LCDM model in the Icr 
CL. 



2. With the current observational data, the variance of wq 
and Wfl we get are still large; the Icr constraints on wq 
and Wa are wo = -1.0537+|J[^32 and w„ = 0.2738+1;;^'^'^. 
This result implies that though the dynamical dark en- 
ergy models are mildly favored, the current data can- 
not distinguish different dark energy models decisively. 
With the fitting results we obtained, we can reconstruct 
the evolution of the EOS of dark energy, w{z) which 
is shown in Fig. [3] From the figure, we can directly 
see that although the quintom model is more favored, 
LCDM, however, still cannot be excluded. 

3. From Fig. |3] one can see that the allowed value of w{z) 



In this paper, we develop a new parameterization which is 
the modification of Efstathiou' parameterization. In this new 
parameterization, there are three free parameters: wq, Wg and 
Q.,n. Then we carry out the global fitting on these model us- 
ing the cun-ent data: SNe la, BAO, CMB and OHD. From the 
analysis, the best fit values of the cosmological parameters 
with Icr confidence-level regions are: Q.,„ - Q.2135^1''l\^^y 
wo = -1.0537+0 [432 and w„ = 0.213?,+lj%l. From the anal- 
ysis, we can directly see that the quintom model is more fa- 
vored, but this result is also consistent with the LCDM model 
in the Icr CL. 
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